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The masses of the low lying baryons are evaluated using two degenerate flavors of twisted mass 
sea quarks corresponding to pseudo scalar masses in the range of about 270 MeV to 500 MeV. The 
strange valence quark mass is tuned to reproduce the mass of the kaon in the physical limit. The 
tree- level Symanzik improved gauge action is employed. We use lattices of spatial size 2.1 fm and 
2.7 fm at two values of the lattice spacing with ro/a = 5.22(2) and ro/a = 6.61(3). We check 
for both finite volume and cut-off effects on the baryon masses. We performed a detailed study of 
the chiral extrapolation of the octet and decuplet masses using SU(2) x?T. The lattice spacings 
determined using the nucleon mass at the physical point are consistent with the values extracted 
using the pion decay constant. We examine the issue of isospin symmetry breaking for the octet and 
decuplet baryons and its dependence on the lattice spacing. We show that in the continuum limit 
isospin breaking is consistent with zero, as expected. The baryon masses that we find after taking 
the continuum limit and extrapolating to the physical limit are in good agreement with experiment. 
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I. INTRODUCTION 

In the last couple of years an intense and successful effort in extending unqucnchcd lattice calculations towards 
realistic values of quark masses, small lattice spacings and large volumes has been undertaken using a variety of 
algorithmic techniques and lattice actions. A review of the salient features of the various discretization schemes 
currently employed can be found in Ref [1]. Of particular relevance to the current work are the calculations of the 
low-lying baryon spectrum using two degenerate flavors (Nf = 2) of light dynamical quarks. Such studies have 
been carried out by the MILC collaboration [2, 3] using Kogut-Susskind fermions and by the European Twisted 
Mass Collaboration (ETMC) [4] for the nucleon (N) and A baryons using twisted mass fermions. There are also 
baryon mass calculations using two degenerate flavors of light quarks and a strange quark with the mass tuned to its 
physical value (Nf = 2+1) mainly using clover improved Wilson fermions with different levels of smearing, such as 
the calculation of the nucleon mass by the QCDSF-UKQCD collaboration [5], and the evaluation of the octet and 
decuplet spectrum by the PACS-CS [6] and BMW [7] collaborations. The LHP Collaboration computed the octet 
and decuplet spectrum using a hybrid action with domain wall valence fermions on asqtad improved staggered sea 
quarks [8] . Preliminary results on the nucleon mass are also computed using Nf = 2 + 1 domain wall fermions by the 
RBC-UKQCD collaboration [9, 10]. 

In this work we study the low-lying spectrum of the baryon octet and decuplet with twisted mass fermions at 
maximal twist. The light quarks arc dynamical degrees of freedom while in the strange sector we use an Osterwalder- 
Seiler valence quark, following the approach employed in the study of the pseudo scalar meson decay constants [11, 12]. 
The bare strange valence quark mass is taken to be the same as the one determined in the meson studies tuned by 
requiring that the mass of the kaon at the physical point matches its physical value. Using the ETMC Nf = 2 
configurations [13, 14] we calculate the baryon spectrum for pion masses in the range of 270 MeV to 500 MeV and 
at two values of the lattice spacing corresponding to (3 = 3.9 and f3 = 4.05 with r^/a = 5.22(2) and r^/a = 6.61(3), 
respectively, where r$ is determined from the force between two static quarks. Results are also obtained at a third 
/3-value, namely (3 = 3.8, which corresponds to r^/a = 4.46(3). The latter results are not taken into account in the 
final analysis due to large autocorrelation effects observed in the Monte Carlo history for quantities like the PCAC 
mass and the plaquette at small sea quark masses. Data at j3 = 3.8 are only used as a consistency check of the 
continuum extrapolation. For the nucleon mass we also performed the calculation at an even finer value of the lattice 
spacing corresponding to r /a = 8.31(5) and (3 = 4.2 to ensure that indeed the continuum extrapolation using a 
weighted average with results at (3 = 3.9 and j3 = 4.05 is valid. We find that the baryon masses considered here show 
a very weak dependence on the lattice spacing and are fully compatible with an 0(a 2 ) behaviour with an almost 
vanishing coefficient of the a 2 term. This justifies neglecting the 0(a 2 ) term in extrapolating results to the continuum 
limit. 

For a fixed value of the lattice spacing we have used up to five different light quark masses and two different volumes. 
The corresponding m^L values are in the range 3.3 to 7.4, where L is the spatial extent of the lattice. Using these 
various values of the lattice spacing, quark masses and volumes allows us to estimate the volume corrections and 
perform a continuum and chiral extrapolation. The good precision of our results on the baryon masses allows us to 
perform a study of chiral extrapolations to the physical point. This study shows that one of our main uncertainties in 
predicting the mass at the physical point is caused by the chiral extrapolations. Another source of systematic error 
is the partially quenched approximation that we have used. 

An important issue is the restoration of the explicitly broken isospin symmetry in the continuum limit. At finite 
lattice spacing, baryon masses display 0(a 2 ) isospin breaking effects. There are, however, theoretical arguments [15] 
and numerical evidences [1, 16] that these isospin breaking effects are particularly pronounced for the neutral pseudo 
scalar mass whereas for other quantities studied so far by ETMC they are compatible with zero. In this paper we 
will demonstrate that also in the baryon sector these isospin breaking effects are in general small or even compatible 
with zero. For a preliminary account of these results see Ref. [17]. 

The paper is organized as follows: The details of our lattice setup, namely those concerning the twisted mass 
action, the parameters of the simulations and the interpolating fields used, are given in Section II. Section III contains 
the numerical results of the baryon masses computed for different lattice volumes, lattice spacings and bare quark 
masses as well as the Gell-Mann Okubo relations that are supposed to be fulfilled in the exact SU(3) limit. Lattice 
artifacts, including finite volume and discretization errors are discussed in Section IV, with special emphasis on the 
0(a 2 ) isospin breaking effects inherent in the twisted mass formulation of lattice QCD. The chiral extrapolations are 
analyzed in Section V. Section VI contains a comparison with other existing calculations and conclusions are finally 
drawn in Section VII. 
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II. LATTICE SETUP 



The lattice action 



For the gauge fields we use the tree- level Symanzik improved gauge action [18], which includes besides the plaquette 
term /T^* 1 ,, also rectangular (1 x 2) Wilson loops U^^ v 



(1) 



l 



1 



with b\ = —1/12 and the (proper) normalization condition bo = 1 ■ 
the usual Wilson plaquette gauge action. 

The fermionic action for two degenerate flavors of quarks in twisted mass QCD is given by 



i. Note that at b\ = this action reduces to 



(2) 



with t 3 the Pauli matrix acting in the isospin space, fi the bare twisted mass and Dy/ the massless Wilson-Dirac 
operator given by 



where 



V M V(z) = - 



D W [U] 



and V*^(a;) 



UUx — ap,)ijj(x — aft) — ip(x) 



(3) 



(4) 



Maximally twisted Wilson quarks are obtained by setting the untwisted quark mass mo to its critical value m cr , while 
the twisted quark mass parameter jj, is kept non-vanishing in order to be away from the chiral limit. In Eq. (2) the 
quark fields \ are m the so-called "twisted basis" . The "physical basis" is obtained for maximal twist by the simple 
transformation 



ip{x) 



CX P I 



1 X(x), 



tp(x) 



In terms of the physical fields the action is given by 



IJ 5 T 



X(x) exp | ^j-75T 



-y V M V* + ro cr ) + n ) ip(x) 



(5) 



(6) 



In this paper, unless otherwise stated, the quark fields will be understood as "physical fields", ip, in particular when 
we define the baryonic interpolating fields. 

A crucial advantage of the twisted mass formulation is the fact that, by tuning the bare untwisted quark mass too to 
its critical value TO cr , physical observables are automatically O(a) improved. In practice, we implement maximal twist 
of Wilson quarks by tuning to zero the bare untwisted current quark mass, commonly called PC AC mass, mpcAc, 
which is proportional to too — to ct up to 0(a) corrections. As detailed in Ref. [19], topcac is conveniently evaluated 
through 



?TlPCAC 



E x (tUj(x,t)P b (o)) 

t/o»i 2J2JP b (x,t)P b (0)) 



lim 



1.2 



(7) 



~, b ~ . b 

where A ^ = XI^I^^X ^ s ^ ne ax i a l vector current and P = X75"Y^ ^ s ^ nc P scu do scalar density in the twisted basis. 
The large t/a limit is required in order to isolate the contribution of the lowest- lying charged pseudo scalar meson 
state in the correlators of Eq. (7). This way of determining topcac is equivalent to imposing on the lattice the validity 
of the axial Ward identity between the vacuum and the charged one-pion zero three-momentum state: 



^ = 2to P cac^ 6 , 6=1,2 



(8) 
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The value of m cr is determined at each (3 value at the lowest twisted mass used in our simulations, a procedure that 
preserves 0(a) improvement and keeps 0{a 2 ) small [14, 20]. The twisted mass fermionic action breaks parity and 
isospin at non- vanishing lattice spacing, as it is apparent from the form of the Wilson term in Eq. (6). In particular, 
the isospin breaking in physical observables is a cut-off effect of 0(a 2 ) [21]. To simulate the strange quark in the 
valence sector several choices are possible. We consider a quenched Ostcrwaldcr-Sciler fermion [22] with the following 
action in the twisted basis: 

S s = a 4, ^2x s {x)(D w [U} + m + i^ s j 5 )xs(x) . (9) 

X 

This is naturally realized in the twisted mass approach by introducing an additional doublet of strange quark and 
keeping only the positive diagonal component of r 3 . The m value is taken to be equal to the critical mass determined 
in the light sector, thus guaranteeing the O(a) improvement in any observable. The reader interested in the advantage 
of this mixed action in the mesonic sector is referred to the Refs [11, 12, 23, 24, 25]. 

B. Simulation details 

The input parameters of the calculation, namely P, L/a and a/i arc summarized in Table I. The corresponding lattice 
spacing a and the pion mass values, spanning a mass range from 270 MeV to 500 MeV, are taken from Ref. [26]. At 
m-a f=a 300 MeV we have simulations for lattices of spatial size L = 2.1 fm and L = 2.7 fm at p = 3.9 allowing to 
investigate finite volume effects. Finite lattice spacing effects are investigated using two sets of results at j3 = 3.9 
and (3 — 4.05. The set at (3 = 3.8 is used only as a cross-check and to estimate cut-off errors. These sets of gauge 
ensembles allow us to estimate all the systematic errors in order to have reliable predictions for the baryon spectrum. 



P = 4.05, a = 0.0666(6) fm from f n [26], r /a = 6.61(3) 


32^ x 64, L = 2.13 fm 


an 0.0030 0.0060 0.0080 0.012 
No. of confs. 269 253 409 182 
m n (GeV) 0.3070(18) 0.4236(18) 0.4884(15) 0.6881(18) 
m^L 3.31 4.57 5.27 7.43 


P = 3.9, a = 0.0855(6) fm, from f w [26], r /a = 5.22(2) 


24 a x 48, L = 2.05 fm 
32 3 x 64, L = 2.74 fm 


an 0.0030 0.0040 0.0064 0.0085 0.010 

No. of confs - 782 545 348 477 

m n (GeV) - 0.3131(16) 0.3903(9) 0.4470(12) 0.4839(12) 

m n L 3.25 4.05 4.63 5.03 

an 0.003 0.004 

No. of confs 659 232 

m n (GeV) 0.2696(9) 0.3082(6) 

m^L 3.74 4.28 


= 3.8, a = 0.0995(7) fm r /a = 4.46(3) 


24 a x 48, L = 2.39 fm 


an 0.0060 0.0080 0.0110 0.0165 
No. of confs 215 302 248 244 

(GeV) 0.3667(17) 0.4128(16) 0.4799(9) 0.5855(10) 
rn^L 4.44 5.00 5.81 7.09 



TABLE I: Input parameters (P,L,n) of our lattice calculation and corresponding lattice spacing (a) and pion mass (m^). 



C. Tuning of the bare strange quark mass 

In a previous paper from the ETM collaboration [11], pseudo scalar meson masses have been computed for different 
values of the sea and valence quark masses for the /3 = 3.9 gauge configurations. Using the experimental value 
of the mass ratio of the kaon to the pion, mK/m^, the bare strange quark mass can be set. We use the value of 
a^s = 0.0217(22) at P = 3.9 taken from Table 2 of Ref. [11]. In a more recent study of the pseudo scalar decay 
constant of kaons and D-mesons [12], the computation was extended to P = 3.8 and P = 4.05. However, this is still a 
preliminary analysis and an ongoing analysis for the accurate extraction of quark masses is still in progress. One can 
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obtain an estimate of the bare strange quark mass at a given value of (3 by taking the results at (3 = 3.9 as a reference 
and using the scaling relation [27] : 



a/j s ((3) 



Z P {d) 



a((3) 



Z p (/3 = 3.9) a(/3 = 3.9) 



a/x s (/3 = 3.9) . 



(10) 



The values we use for (3 = 3.8 and (3 = 4.05 given in Table II are obtained by applying Eq. (10). We use the value of 
the renormalization constant Z p {(3) found in the preliminary analysis of Ref. [28] within the RI'-MOM scheme. This 
value is in agreement with a complementary analysis given in Ref. [29]. 



(3 = 3.8 (3 = 3.9 (3 = 4.05 

0.0208(15) (48) 0.0217(22) 0.0166(18)(29) 



TABLE II: Bare strange quark mass used in the valence sector for different (3 values. 



D. Interpolating fields 

The low lying baryons belonging to the octet and decuplet representations of SU(3) are given in Figs. 1 and 2 
respectively. They are classified by giving the isospin, /, the third component of the isospin, 1^, the strangeness 
(s), spin and parity. In order to extract their masses in lattice QCD we evaluate two point correlators. We use 
interpolating fields to create these states from the vacuum that have the correct quantum numbers and reduce to the 
quark model wave functions in the non-relativistic limit. The interpolating fields used in this work are collected in 
Tables III [30, 31] and IV [30, 32] for the octet and decuplet respectively. 



s = 



s = -1 £- 



s = -2 




A 




E+ 



q = l 



5 = 



S= -1 



s = -2 



s = -3 



A+ A+^ 




q = 2 



q=-l 



q = -l 



FIG. 1: The low lying baryons belonging to the octet rep- 
resentation labeled by the value of ^3 and hypercharge. 



FIG. 2: The low lying baryons belonging to the decu- 
plet representation labeled by the value of I3 and hyper- 
charge. 



Local interpolating fields are not optimal for suppressing excited state contributions. We instead apply Gaussian 
smearing to each quark field, g(x, t): g smear (x, t) = J2 y ^( x ! y; U (t))q(y, t) using the gauge invariant smearing function 



F(x,y;U(t)) = (l+aH) n (x >T) U(t)), 



constructed from the hopping matrix, 



iJ(x, y; U(t)) = J2 (W, *)<W-i + # (x - i, t)-5x, y+J ) ■ 
i=l ^ ' 



(11) 



(12) 



Furthermore we apply APE smearing to the spatial links that enter the hopping matrix. The parameters of the 
Gaussian and APE smearing are the same as those used in our previous work devoted to the nucleon and A masses [4] . 
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Strangeness 


Baryon 


Interpolating field 


I h 


s = 


P 




1/2 +1/2 


n 


X n = ^abc(dlC-y 5 Ub)d c 


1/2 -1/2 




A 


X* 8 = -^tabc{'2{ulC>ysd b )s c + (ulC'f5Sb)d c - (dlC"/ 5 s b )u c } 





s = 1 


E+ 


X S = £abc{u r C'y5Sb)u c 


1 +1 




E° 


X S = -^^abc{(u!C'yaSb)dc + (daC^f5Sb)Uc} 


1 +0 




£- 


X S = eabc(d% C"/ 5 Sb)d c 


1 -1 


s = 2 




X" ~ £abc{S a Cy5Ub)Sc 


1/2 +1/2 






X~ = £abc(saC-y 5 d b )s c 


1/2 -1/2 



TABLE III: Interpolating fields and quantum numbers for the baryons in the octet representation. 



Strangeness 


Baryon 




Interpolating field 


I h 




A++ 




A++ / Try \ 
Xm = ^abc{U a CllJ.Ub)Uc 


3/2 +3/2 


s = 


A+ 
A 




■^^abc{2(UaC'y^db)u c + («I C7 M u 6 )d c } 
^=e a6c {2(da C7 M iti,)d c + {dlC^^d b )u c } 


3/2 +1/2 
3/2 -1/2 




A" 




Xm = £abc(daC"/ ti db)d c 


3/2 -3/2 




E*+ 


£*+ _ 


-^^abc{(UaC-y^Ub)s c + 2(s„ Cj^Ub )u c } 


1 +1 


s = 1 


E*° 


Aji ^ V 3 t abc 


{(«! Cj^d b )s c + (daCj^Sb)^ + (s^ C7 M M6)d c } 


1 +0 




E*~ 


X(U 


^ea6c{(daC7 fl di,)s c + 2(s^C7 M d 6 )d c } 


1 -1 


s = 2 


=*o 




Xm ° = e a(, c (sa C%Uj,)s c 


1/2 +1/2 








Xm = eab c (srC7 M di,)s c 


1/2 -1/2 


s = 3 






Xm = tabcislC-y^s^Sc 


+0 



TABLE IV: Interpolating fields and quantum numbers for baryons in the decuplet representation. 



E. Two-point correlators 



To extract masses in the rest frame we consider two-point correlators defined by 

Cx(t>P = 0) = = ' = ^ 4 ) (■/x( x sink, isink)-7x( x sourcej ^source)}) t — t s ink — ^source • (13) 

x aink 

Space-time reflection symmetries of the action and the anti-periodic boundary conditions in the temporal direction 
for the quark fields imply, for zero three- momentum correlators, that Ci(t) = —C^(T — t). So, In order to decrease 
errors we average correlators in the forward and backward direction and define: 

C x (t) = C x (t)-C x (T-t). (14) 

In order to decrease correlation between measurements, we choose the source location randomly on the whole lattice 
for each configuration. Masses are extracted from the so called effective mass which is defined by 

/ l + V°° c e Ait \ 

am^(t) = - \og{C x (t)/C x (t - 1)) = am x + log - ^j£ i=1 * _„ — > am x , (15) 

where A j = — m x is the mass difference of the excited state i with respect to the ground state mass m x . 

In Figs. 3 and 4 we show the effective masses of the baryons in the octet and decuplet representation respectively. 
As can be seen a plateau region can be identified for all baryons. What is shown in these figures are effective masses 
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extracted from correlators where smearing is applied both at the sink and source. Although local correlators are 
expected to have the same value in the large time limit, smearing suppresses excited state contributions yielding a 
plateau at earlier time separations and a better accuracy in the mass extraction. Our fitting procedure to extract mx 
is as follows: The mass is obtained from the leading term in Eq. (15), i.e. from a constant fit to mx- A second fit, 
including the first excited state, allows us to estimate the systematical error of the previously determined mx due to 
excited states for a given plateau range. The plateau range is then chosen such that the systematical error on mx 
drops below 50% of its statistical error This criterion is in most of the cases in agreement with a x 2 /d.o.i. < 1. In 
the cases in which this criterion is not satisfied a careful examination of the effective mass is made to ensure that the 
fit range is in the plateau region. The results for the masses of the octet and decuplet at j3 = 3.9 are are collected in 
Tables V and VI respectively. The corresponding results for the masses at j3 = 4.05 are given in Table VII and VIII. 
The errors are evaluated using both jackknife and the T-method [33] to ensure consistency. 



3 I 
Ai 
Nl 



I I 
I 



1 ; 1 1 1 

1 1 1 1 1 1 1 



t/a 



FIG. 3: Effective masses of the octet states for /3 = 3.9, 
afi = 0.004 on a 32 3 x 64 lattice using 232 configurations. 



0.6 



a 
z 

A 



A 



t/a 



FIG. 4: Effective masses of the decuplet states for (3 = 
3.9, afi = 0.004 on a 32 3 x 64 lattice using 232 configu- 
rations. 



III. RESULTS 



The bulk of the numerical results are presented in this section. Baryon masses are given in lattice units. Our 
procedure to convert the results to physical units will be discussed in the next section. 



A. Baryon masses 



In Tables V to VIII we present the masses of the octet and decuplet states with the lattice input parameters given in 
Table I. For the isospin multiplcts we have computed separately the masses corresponding to each isospin components 
as well as their averaged value. These results (averaged values in case of isospin multiplcts) are displayed in Figs. 5 
and 6. The j3 = 3.9, L/a = 24 data are linked by dotted lines to guide the eye. An inspection of the plots indicates 
that the lattice artifacts, studied in detail in the next section, are small. Notice that the natural order of the S* and 
S* states comes out to be correct for m v < 300 MeV while for larger masses this order is inverted. 
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Ufl 


stat. 


amjv omA am E Av am E + sm E « am E - am 3 A, am H o am H _ 


24 a x 48 


0.0040 
0.0064 
0.0085 
0.0100 


782 
545 
348 
477 


0.5111(58) 0.5787(42) 0.6075(46) 0.6175(66) 0.6118(48) 0.5959(52) 0.6497(31) 0.6695(42) 0.6372(31) 
0.5514(49) 0.6017(42) 0.6265(48) 0.6487(72) 0.6278(52) 0.6131(52) 0.6636(36) 0.6876(50) 0.6500(36) 
0.5786(67) 0.6198(51) 0.6491(55) 0.6679(62) 0.6529(49) 0.6358(46) 0.6728(43) 0.6956(58) 0.6593(43) 
0.5973(43) 0.6326(36) 0.6522(41) 0.6662(56) 0.6539(43) 0.6429(44) 0.6793(36) 0.6959(49) 0.6683(32) 


32 a x 64 


0.0030 
0.0040 


652 
232 


0.4958(43) 0.5613(33) 0.5891(42) 0.6069(68) 0.5932(50) 0.5775(39) 0.6382(30) 0.6572(44) 0.6275(26) 
0.5126(46) 0.5750(35) 0.6117(40) 0.6281(73) 0.6158(40) 0.5960(47) 0.6511(34) 0.6748(46) 0.6358(32) 



TABLE V: Baryon masses in the octet representation at /3 = 3.9 in lattice units. 



Ufl 


stat. 


am A ++,~ 


am A +,o 


om E »Av am E * + 


am s .o 


am E .- 


am H «Av am s .o am s .- omn 


24 d x 48 


0.0040 
0.0064 
0.0085 
0.0100 


782 
545 
348 
477 


0.660(14) 
0.709(11) 
0.714(12) 
0.7531(67) 


0.670(13) 
0.711(12) 
0.733(13) 
0.7559(75) 


0.7166(82) 0.709(11) 
0.7461(84) 0.740(10) 
0.7517(88) 0.739(11) 
0.7794(66) 0.7808(62) 


0.7226(81) 
0.7480(93) 
0.760(11) 
0.7809(64) 


0.7222(95) 
0.7489(93) 
0.7645(98) 
0.7798(69) 


0.7311(51) 0.7381(59) 0.7200(66) 0.8079(52) 
0.7412(78) 0.7552(76) 0.7344(84) 0.8156(63) 
0.7415(85) 0.7529(81) 0.7367(84) 0.8133(66) 
0.7618(73) 0.7741(64) 0.7484(74) 0.8284(51) 


Z2 6 x 64 


0.003 
0.004 


652 
232 


0.6234(139) 0.6497(133) 0.6859(96) 0.6838(93) 0.6859(106) 0.7027(101) 0.7058(50) 0.7097(58) 0.7032(53) 0.7926(49) 
0.651(16) 0.659(15) 0.713(10) 0.705(12) 0.716(12) 0.7173(99) 0.7291(74) 0.7366(79) 0.7192(72) 0.8037(69) 



TABLE VI: Baryon masses in the decuplet representation at /3 = 3.9 in lattice units. 



a\i stat. arriN amA amsav am s + am E « am s - am=a^ am 3 o am s - 

32 a x 64 

0.0030 269 0.4091(60) 0.4540(38) 0.4761(44) 0.4885(62) 0.4774(47) 0.4651(53) 0.5082(31) 0.5177(39) 0.5007(29) 

0.0060 253 0.4444(47) 0.4792(47) 0.4944(44) 0.5022(66) 0.4960(45) 0.4834(45) 0.5192(42) 0.5277(50) 0.5112(37) 

0.0080 409 0.4714(31) 0.4957(30) 0.5089(31) 0.5179(41) 0.5095(32) 0.5019(31) 0.5262(28) 0.5350(34) 0.5199(25) 

TABLE VII: Baryon masses in the octet p representation at (3 = 4.05 in lattice units. 





stat. 


am A ++,- am A +,o ams-av am E *+ am s *o am E »- ams*av am=.o am 3 »- amn 


32 J x 64 


0.0030 
0.0060 
0.0080 


269 
253 
409 


0.5381(93) 0.5441(93) 0.5728(79) 0.5673(94) 0.5750(86) 0.5734(80) 0.5772(56) 0.5796(54) 0.5750(56) 0.6361(46) 
0.5505(77) 0.5581(90) 0.5805(66) 0.5754(71) 0.581(11) 0.5844(68) 0.5816(47) 0.5834(50) 0.5802(48) 0.6286(53) 
0.5918(60) 0.5906(63) 0.6078(59) 0.6044(68) 0.5850(74) 0.6099(57) 0.5940(43) 0.6021(50) 0.5873(43) 0.6461(49) 



TABLE VIII: Baryon masses in the decuplet representation at /3 = 4.05 in lattice units. 



B. Strange quark mass dependence 

The dependence of the masses of baryons with strangeness on the bare strange quark mass has been investigated 
at (3 = 3.9 for a/j, = 0.004. The results are given in Tables IX and X and displayed in Figs. 7 and 8. The vertical 
dotted line indicates the value of the tuned bare strange quark mass as given in Table II. The SU (3) symmetric point 
Us = \i is given by the nucleon and A mass for the octet and decuplet respectively. As can be seen in the SU (3) 
limit all the octet and decuplet masses converge to a single point up to cut-off effects and the fact that we only have 
Nf = 2 simulations. For clarity we only show in Fig. 7 the mass of A, S Av and S Av . They should be degenerate with 
the nucleon in the limit of fi s = fi. Indeed, if one computes the nucleon mass with the same statistics with that used 
for E Av and S Av , one finds them to be degenerate within the errors as can be seen in Fig. 7. 

The corresponding results for the decuplet-baryons are displayed in Fig. 8. As can be seen, also in the case of the 



9 



I 
A 



> (J=3.9 L=24 : 
P=3.9 L=32 = 
P=4.05 L=32 i 



0.3 

iti (GeV) 



FIG. 5: Octet states measured in our different gauge en- 
sembles. Physical points are indicated by their name 
symbols (magenta). The data at /3 = 3.9, L/a — 24 are 
connected by dotted lines to guide the eye. 
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FIG. 6: The same as Fig. 5 but for the decuplet states. 



decuplet masses there is convergence to the A mass as predicted in the exact SU (3) limit [i s = /i. 




0.05 



FIG. 7: Masses for octet baryons at j3 = 3.9 and a/x = 
0.004 on a lattice of size 24 3 x 48 as a function of afi B . 
The vertical dashed line indicates the value of the tuned 
bare strange quark mass. The dotted lines are to guide 
the eye. 
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FIG. 8: The same as for Fig. 7 but for the decuplet 
baryons. 



The ii s dependence of the strange baryon masses provides an estimate of systematic errors due to the uncertainty 
in the tuning of the strange quark mass. As already explained, the kaon mass at the physical point is used to fix /j, s . 
This gives a/j, s = 0.0217(22). The ~ 10% uncertainty leads to a corresponding error in the strange baryon masses 
that can be estimated by the variation of their masses in the vicinity of fj, s . At \i = 0.004 we estimate an error that 
is comparable to the statistical error. In what follows we will analyze our results taking into account only statistical 
errors. This analysis shows that when the statistical error is given on the final results of strangeness non-zero baryon 
masses one must bear in mind that there is a systematic error of about the same magnitude due to the strange quark 
mass determination. 
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stat. 


amA 


am E Av 


am E + am E o 




am g i. 


am H o 


am 3 _ 


24 3 x 48 



0.0064 597 0.533(8) 0.545(7) 0.549(12) 0.563(5) 0.537(6) 0.545(7) 0.560(11) 0.530(6) 

0.0085 316 0.537(10) 0.557(11) 0.559(19) 0.557(12) 0.554(8) 0.563(9) 0.585(14) 0.549(8) 

0.0100 316 0.542(9) 0.564(10) 0.567(18) 0.564(11) 0.561(7) 0.574(8) 0.597(13) 0.560(7) 

0.0175 315 0.563(8) 0.596(9) 0.600(14) 0.593(9) 0.593(7) 0.626(6) 0.644(8) 0.610(6) 

0.0200 308 0.568(8) 0.606(8) 0.609(13) 0.602(9) 0.603(6) 0.641(6) 0.660(8) 0.625(5) 

0.0250 311 0.584(7) 0.626(6) 0.627(12) 0.619(8) 0.620(6) 0.671(5) 0.688(7) 0.656(5) 

0.0400 316 0.624(7) 0.674(6) 0.676(10) 0.667(7) 0.672(6) 0.751(4) 0.764(5) 0.738(4) 

0.0800 314 0.718(7) 0.780(5) 0.787(7) 0.772(7) 0.776(7) 0.935(3) 0.945(4) 0.926(3) 

TABLE IX: Octet masses for (3 = 3.9, a/i = 0.004 on a 24 3 x 48 lattice as a function of a/i s . 



afi a stat. ojn s ,Av am E * + am E .o am E »- nm 3 ,Av am H *o am s ._ amn 

24 a x 48 

0.0064 597 0.665(12) 0.658(18) 0.669(14) 0.669(14) 0.636(9) 0.645(12) 0.628(9) 0.678(14) 

0.0085 316 0.695(19) 0.719(13) 0.713(16) 0.733(9) 0.648(11) 0.670(15) 0.624(10) 0.734(11) 

0.0100 316 0.700(18) 0.722(12) 0.715(15) 0.697(22) 0.658(9) 0.680(13) 0.636(9) 0.744(10) 

0.0175 315 0.721(14) 0.729(14) 0.736(12) 0.718(19) 0.705(7) 0.722(8) 0.690(7) 0.796(6) 

0.0200 308 0.725(14) 0.734(13) 0.741(11) 0.718(17) 0.720(6) 0.734(7) 0.704(7) 0.807(7) 

0.0250 311 0.740(13) 0.735(16) 0.753(11) 0.740(17) 0.747(6) 0.759(6) 0.733(6) 0.838(6) 

0.0400 316 0.778(11) 0.770(13) 0.788(9) 0.778(15) 0.821(5) 0.831(5) 0.811(5) 0.934(4) 

0.0800 314 0.864(7) 0.854(11) 0.870(8) 0.865(9) 0.993(4) 0.996(5) 0.987(4) 1.169(3) 

TABLE X: Decuplet masses for (5 = 3.9, cifi = 0.004 on a 24 3 x 48 lattice as a function of afi s 



C. Gell-Mann-Okubo relation 



Assuming a small SU(3) breaking, Okubo derived interesting relations among baryons masses. We examine in this 
section how well the Gcll-Mann-Okubo (GMO) relations [34] are fulfilled for the baryons masses obtained on our 
lattices at different pion mass values. As we will discuss in detail in the next Section, volume and discretization effects 
are small, and therefore it suffices to analyze the (3 = 3.9 and L = 24 x 48 results. For this study we use the lattice 
spacing determined from to convert to physical units. 

For the J p = l/2 + octet the GMO relation can be written in the form: 

Mh + Mjv 3M a + M e 

2 = 4 ' (16) 

The results are displayed in Fig. 9 where the left and right hand side terms of Eq. (16) are separately plotted as a 
function of mj. The difference between the two terms are compatible with zero at any pion mass. The experimental 
values, shown by the squares, are respectively 254 MeV and 248 MeV. These results are similar to those presented in 
Ref. [35] using a mixed action setup with valence domain wall fermions on rooted staggered sea fermions. 

For the J p = 3/2+ decuplet, the GMO relations predict equal mass difference among two consecutive (AS = 1) 
isospin multiplets: 

M s . - A/ A = Ms* - M E . = M Q - M E . . (17) 

The results for the decuplet baryons are displayed in Fig. 9. As can be seen, the equalities of Eq. (17) are strongly 
violated; the three mass differences of Eq. (17) are spread over about 200 MeV for the range of pion masses that have 
been computed. The experimental values for these mass differences are 153, 149, 139 MeV, shown in the plot by the 
squares. In the lattice results the larger deviation comes from Ms* — Ms* , while for M E * — Ma and Mq — M=» the 
mass differences are smaller. The mass difference Af~» — A/s* is increasing as the pion mass decreases. Unfortunately, 
with our present statistics it is unclear whether this increase is sufficient to bring this mass difference in agreement 
with experiment but the trend is definitely in the right direction. 

A third relation exists, that connects the J p = 1/2+ octet masses with the J p — 3/2 + decuplet masses, which 
reads as: 



3M A - M s - 2M N = 2(M S . - M A ) 



(18) 
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Experimentally, this relation is fulfilled at the 10% level yielding 276 MeV for the left hand side and 305 MeV for the 
right hand side of Eq. (18). These values are again shown by the filled squares in Fig. 9. The corresponding lattice 
results are shown in the same figure. One can see that, as in the octet case, the relation of Eq. (18) is satisfied within 
our statistical uncertainties at each pion mass. It also approaches the experimental results with decreasing pion mass. 
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FIG. 9: Gell Mann Okubo relations for the baryon octet (a) decuplet (b) and mixed octet-decuplet (c) as a function of m\. 
Vertical lines correspond to the physical results. Data were obtained from simulations at /3 = 3.9 and L/a — 24. The lattice 
spacing determined from is used to convert to physical units. 



Fulfillment of the GMO relations is considered a success of SU(3) symmetry. Violations of these relations indicate 
that SU (3) breaking is not small. Therefore one would expect that these relations are better satisfied as we approach 
the SU(3) limit a[i = afi s = 0.0217, up to discretization effects. This corresponds to about m 2 ~ 0.50 GeV 2 . For the 
decuplet mass relation given in Eq. (17) it is unclear if this would be indeed satisfied by the lattice data whereas the 
other two relations arc fulfilled at all masses. 



IV. SYSTEMATICS 



In order to compare our lattice results collected in Tables V, VI, VII and VIII to the physical masses we need to 
check for finite volume effects, cut-off effects and the extrapolation to the physical light quark masses. The strange 
quark was fixed to the physical value using the kaon mass with the light quarks extrapolated to the physical point as 
explained in Section II. B. A check of the effect of this tuning on baryon masses has been discussed in Section III.B. 
In this section we discuss finite volume and cutoff-effects, in particular the isospin breaking. 



A. Finite volume effects 



Finite volume corrections to the nucleon mass in Nf = 2 lattice QCD have been studied in Ref. [5] within the p 
expansion which assumes that finite size effects originate from pions that propagate around the spatial box. Using 
relativistic SU(2) baryon chiral perturbation theory [36] the finite volume corrections to the nucleon mass to 0(p 4 ) 
are: 



mjv(oc) = itin(L) — 8m a (L) — drrib(L) 



(19) 
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where 



Sm a (L) 
6m b (L) 



3g 2 A m° N ml 



3m 



8n 2 f, 



2 l-co ' / 



f2 



(2ci 



C3 



K 1 {L\n\m v ) K 2 (L\n\m w ) 



L|n|m„ 



C2- 



(L|n|m„ 



(20) 



K v (x) is the modified Bessel function and the sum is over all integer vectors n excluding n = 0. The parameters vnP N 
and ci are determined by fitting first the nucleon mass to the same order [37, 38, 39] given by 



m N 



3fl 



167T/2 



4Si(A)m* 



37«: 



16^ 2 /2 



f'2 



25| 
o 

AT 



III 



c 2 — 8ci + 4c 3 



'■jV 



log 



A 



(21) 



We take the cut-off scale A = 1 GcV, f„ = 130.70 MeV and fix the dimension two low energy constants c 2 = 
3.2 GeV -1 [40] and C3 = —3.45 GcV" 1 [38, 41]. These values are consistent with empirical nucleon-nucleon phase 
shifts [42, 43]. The counter-term E\ is taken as an additional fit parameter. We then use these parameters to estimate 
the volume corrections to the nucleon mass. The results are listed in Table XI. As can be seen the corrections for our 
lattices are, in all cases except one, smaller than the statistical errors. In the analysis that follows we will use the 
volume corrected nucleon mass. 



am-r: 


amjv(L) aS a (L) + a5b(L) 




= 3.9 24 a 


x 48 


0.1362 


0.5111(58) 


0.0068 


0.1684 


0.5514(49) 


0.0046 


0.1940 


0.5786(67) 


0.0026 


0.2100 


0.5973(43) 


0.0021 




P = 3.9 32 3 


x 64 


0.1168 


0.4958(34) 


0.0014 


0.1338 


0.5126(46) 


0.0011 




f3 = 4.05 32 3 


x 64 


0.1038 


0.4091(60) 


0.0035 


0.1432 


0.4444(47) 


0.0018 


0.1651 


0.4714(31) 


0.0012 



TABLE XI: Volume correction to the nucleon mass. 



Concerning the other baryons a recent analysis using SU(3) heavy baryon chiral perturbation theory has shown that 
the volume corrections are smaller than for the nucleon [44] . Given that the volume correction found for the nucleon 
are smaller than the statistical errors we can safely neglect any volume corrections for the other baryons computed in 
this work. This is also corroborated by our lattice results at a,[i = 0.004 where simulations at two volumes are used. 



B. Isospin breaking 

The twisted mass action breaks isospin explicitly to 0(a 2 ). How large this breaking is depends on the size of the 
0(a 2 ) terms. It was shown that this cut-off effect is particularly large for the neutral pion [15] but small for other 
quantities. Indeed we verified that isospin breaking between the A ++, ° and A +,_ is consistent with zero for lattice 
spacings below about 0.1 fm [4]. We here address this issue for the octet and decuplet baryons. We show in Fig. 10 
the mass differences for the E, S, A, E* and S* charge multiplets as a function of the pion mass at two values of /3. 

As can be seen, we confirm that for the A system isospin breaking is consistent with zero. This is also true for the 
E* and S* as well as for the E at the smaller lattice spacing. In the case of the S we observe a non-zero splitting 
that decreases with the lattice spacing. If one interpolates the results at the two (3 values to the same pion mass, as 
discussed in more detail in the next Section, and makes a linear extrapolation in a 2 one finds that this splitting goes 
to zero in the continuum limit as expected. Whereas this confirms that this splitting is a cut-off effect, to perform a 
proper analysis one would need results at an additional lattice spacing. For the current work we conclude that isospin 
splitting at these two lattice spacings is negligible for all baryons expect for the 3 where an isospin breaking of about 
6% is observed that vanishes at a rate proportional to a 2 . 
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C. Continuum extrapolation 



In order to assess cut-off effects we use results at (3 = 3.9 and (3 = 4.05. The lattice results, expressed in units of 
the Sommer scale r . are interpolated to the same pion mass in units of r at each f3- value. We give the interpolated 
results at six values of m^ro in Tables XII and XIII. Interpolating linearly or with one-loop order chiral perturbation 
theory gives values consistent within error bars. Given the size of these cut-off effects a weighted average of the baryon 
masses between these two (3 values gives an estimate of the values in the continuum limit. It must be stressed that 
estimating the strange quark mass at (3 ~ 4.05 using Eq. (10) may cause residual cut-off effects on the few percentage 
level that are not taken into account with the continuum extrapolation as performed here. The results obtained from 
the weighted averaging of data at (3 = 3.9 and (3 = 4.05 are listed in Tables XII and XIII and are plotted in Figs. 11 
and 12 . In the figures we also include results at (3 = 3.8. If cut-off effects are small for all (3- values then results at 
(3 = 3.8 should fall onto the same line. As can be seen this is best fulfilled for the nucleon mass. Furthermore for the 
nucleon and the A we also show results at a smaller value of the lattice spacing corresponding to (3 = 4.2. Essentially, 
the a 2 dependence of the nucleon and A mass as computed at (3 = 3.9, 4.05 and 4.2 is consistent with a constant 
behaviour, verifying that for lattice spacings below 0.1 fm cut-off effects are indeed small. For the A mass results at 
(3 = 3.8, 3.9 and [3 = 4.05 are consistent with a constant. This holds approximately also for the other baryons. Within 
the statistical errors one therefore concludes that for lattice spacings below 0.1 fm cut-off effects are under control. 
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romjv 



r mA 



rom E A» 



rom s + rom E o rom E _ rom H Av. 



rom s o 



rom = ~ 



0.60 


2.571(23) 


2.922(18) 


3.062(22) 


u. / u 


9 £71 (OA\ 


nm ( i 

O.UUl ^lo ) 




0.80 


2.757(30) 


3.085(25) 


3.212(29) 


0.90 


2.880(26) 


3.156(22) 


3.286(25) 


1.00 


2.992(35) 


3.226(27) 


3.381(29) 


1.10 


3.111(23) 


3.305(190 


3.405(21) 


0.60 


2.600(43) 


2.946(29) 


3.107(32) 


0.70 


2.692(40) 


3.008(25) 


3.152(29) 


0.80 


2.788(46) 


3.074(32) 


3.200(34) 


0.90 


2.874(32) 


3.135(32) 


3.242(30) 


1.00 


2.984(33) 


3.205(32) 


3.298(30) 


1.10 


3.119(21) 


3.283(20) 


3.370(21) 


0.60 


2.577(20) 


2.929(15) 


3.077(18) 


0.70 


2.676(21) 


3.003(15) 


3.179(17) 


0.80 


2.766(25) 


3.080(20) 


3.207(22) 


0.90 


2.878(20) 


3.149(18) 


3.267(19) 


0.10 


2.988(24) 


3.217(21) 


3.342(21) 


1.10 


3.116(15) 


3.295(14) 


3.387(15) 



P = 3.9 
3.156(36) 3.084(26 
3.279(38) 3.215(21 
3.318(42) 3.219(31 
3.404(38) 3.293(27 
3.482(33) 3.401(26 
3.477(29) 3.414(22 

(3 = 4.05 
3.199(46) 3.115(34 



3.233(41) 
3.269(49) 
3.295(45) 
3.348(45) 
3.430(27) 



3.161(31 
3.209(36 
3.253(31 
3.308(31 
3.373(21 



continuum limit 
3.173(28) 3.095(21 
3.258(28) 3.198(17 
3.297(32) 3.215(23 
3.359(29) 3.275(20 
3.436(26) 3.363(20 
3.452(20) 3.392(15 



3.004(21 
3.111(25 
3.144(31 
3.215(27 
3.310(24 
3.358(23 

3.034(38 
3.080(35 
3.127(40 
3.165(31 
3.230(31 
3.325(21 

3.011(18 
3.101(20 
3.138(24 
3.193(20 
3.280(19 
3.340(15 



3.324(16 
3.399(18 
3.433(21 
3.472(19 
3.507(23 
3.547(19 

3.336(23 
3.362(21 
3.391(25 
3.418(28 
3.448(29 
3.481(19 

3.328(13 
3.383(13 
3.415(16 
3.456(16 
3.484(18 
3.514(13 



3.421(23 
3.523(24 
3.555(30 
3.598(26 
3.629)30 
3.633(26 

3.400(29 
3.425(26 
3.451(32 
3.474(34 
3.504(34 
3.539(23 

3.413(18 
3.477(18 
3.506(22 
3.552(21 
3.573(23 
3.580(17 



3.271(14) 
3.319(17) 
3.363(21) 
3.401(19) 
3.436(23) 
3.491(17) 

3.287(22) 
3.313(19) 
3.340(24) 
3.364(25) 
3.397(25) 
3.440(17) 

3.275(12) 
3.316(13) 
3.353(16) 
3.387(15) 
3.418(17) 
3.465(12) 



TABLE XII: Octet masses computed at reference pion masses in units of ro and the corresponding continuum limit values. 



r m A Av 



r m E »Av. r m E ,+ r m E *o r m E * 



rom H ,Av 



rom 3 ,d rom H . 



r mn 



0.60 


3.312(51) 


3.565(51) 


3.558(49) 


0.70 


3.419(57) 


3.719(55) 


3.679(65) 


0.80 


3.631(48) 


3.850(50) 


3.824(62) 


0.90 


3.727(42) 


3.907(44) 


3.871(55) 


1.00 


3.761(47) 


3.911(46) 


3.839(55) 


1.10 


3.950(27) 


4.075(35) 


4.085(33) 


0.60 


3.524(47) 


3.769(57) 


3.732(67) 


0.70 


3.581(43) 


3.789(52) 


3.752(62) 


0.80 


3.615(50) 


3.808(60) 


3.773(71) 


0.90 


3.617(40) 


3.804(45) 


3.768(49) 


1.00 


3.718(41) 


3.876(46) 


3.844(49) 


1.10 


3.922(40) 


4.028(39) 


4.007(45) 



0.60 
0.70 
0.80 
0.90 
1.00 
1.10 



(3 = 3.9 
3.564(56) 3.660(53) 
3.735(61) 3.744(52) 
3.856(55) 3.853(55) 
3.918(49) 3.924(49) 
3.954(58) 3.981(52) 
4.081(34) 4.074(36) 

(5 = 4.05 



3.787(65) 
3.803(57) 
3.819(70) 
3.834(76) 
3.851(76) 
3.868(49) 



3.766(58) 
3.795(53) 
3.822(61) 
3.829(46) 
3.901(47) 
4.042(38) 



continuum limit 

rorriir r m A+ +, r m A +, r m E »Av. r m E .+ r m E .o r m E . 



3.671(27) 
3.805(39) 
3.856(46) 
3.873(41) 
3.862(45) 
3.981(38) 

3.806(40) 
3.817(37) 
3.828(43) 
3.829(32) 
3.862(32) 
3.931(29) 

rom H ,Av. 



3.689(31) 
3.845(41) 
3.925(45) 
3.947(40) 
3.922(43) 
4.046(34) 

3.823(39) 
3.832(36) 
3.842(42) 
3.834(34) 
3.882(35) 
3.987(33) 



3.662(28) 
3.754(38) 
3.812(49) 
3.839(44) 
3.840(44) 
3.909(39) 

3.789(40) 
3.802(37) 
3.816(43) 
3.825(33) 
3.847(33) 
3.885(29) 



rom 2 «o rom s . 



4.131(26) 
4.195(36) 
4.252(37) 
4.259(33) 
4.240(35) 
4.328(27) 

4.221(34) 
4.202(30) 
4.183(37) 
4.141(36) 
4.171(36) 
4.277(33) 



3.439(35) 
3.520(49) 
3.623(35) 
3.668(29) 
3.735(30) 
3.935(20) 



3.656(38) 
3.755(38) 
3.833(38) 
3.856(32) 
3.893(32) 
4.054(26) 



3.619(40) 
3.717(45) 
3.802(47) 
3.813(36) 
3.842(37) 
4.058(27) 



3.660(42) 
3.771(42) 
3.841(43) 
3.893(41) 
3.916(46) 
4.013(28) 



3.709(39) 
3.768(37) 
3.839(41) 
3.874(34) 
3.937(35) 
4.059(26) 



3.712(22) 
3.811(27) 
3.841(31) 
3.845(25) 
3.862(26) 
3.949(23) 



3.741(24) 
3.838(27) 
3.880(30) 
3.882(26) 
3.898(27) 
4.016(24) 



3.704(23) 
3.778(26) 
3.814(32) 
3.830(26) 
3.845(26) 
3.893(23) 



4.164(21) 
4.199(23) 
4.217(26) 
4.205(24) 
4.207(25) 
4.307(21) 



TABLE XIII: The same as Table XII but for the decuplet baryons. 



D. Fixing the lattice spacing 



In order to convert to physical units we need to fix the lattice spacing. The values of the lattice spacing given in 
Table I were extracted using the pion decay constant. Equivalently one can determine the value of r by extrapolating 
the results to the physical point. The value obtained is ro = 0.439(25) fm determined in the light meson sector [45] 
where the systematic error is added to the statistical one. Knowing ro and the ratio r^/a one can determine the 
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FIG. 11: Constant extrapolation to the continuum limit 
for the octet baryons. Stars are for rom^ = 0.615, filled 
triangles for rom, = 0.7, open circles for rom T = 0.8, 
open triangles for rom T = 0.9, filled circles for rom^ — 
1.0, rhombii for r^m-n = 1.1 and crosses for rom^ = 1.25. 
The open squares show the extracted continuum value. 
For the nucleon we also show results at j3 = 4.2. For the 
E+ and H° we omit the case rom^ = 1.0 and shift results 
at rom, = 1.25 and rom^ = 0.8 for clarity. 
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FIG. 12: Constant extrapolation to the continuum for 
the decuplet. The notation is the same as in Fig. 11. 
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FIG. 13: Determination of rg with a simultaneous fit to the lattice data at j3 = 3.9 and j3 = 4.05. The asterisks denote the 
physical point determined by the value of by using 0(p 3 ) and C(p 4 ) x?T as described in Ref. [4]. 

The nucleon mass can be used to set the scale and this determination seems natural if one is interested in the study 
of the baryon spectrum. Our data at j3 = 3.9 and (3 = 4.05 do not show significant lattice spacing effects. Therefore we 
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can make a combined fit using continuum chiral perturbation theory results to determine the j-jj . Chiral corrections to 
the nucleon mass are known to 0(p 4 ) within several expansion schemes. We use the same schemes used in Ref. [4, 46]. 
The fits are shown in Fig. 13. Using the 0(p 3 ) result, which is well established, to extrapolate to the physical pion 
mass we obtain Tq = 0.465(6) fm. If we instead use the data at j3 = 3.9 and [3 = 4.05 to perform the continuum 
limit as discussed in the previous subsection and then fit, we find Tq = 0.471 ± Q.006(stat.) ± 0.015(sj/st.) fm. The 
systematic error is due to the interpolation to a fixed value of m^ro and it is obtained by comparing the value of ?*o 
obtained when linear interpolation is used to the one obtained using 0(p 3 ). Furthermore, we take the difference in the 
value of ro obtained using continuous results and the value found by fitting the lattice data at (3 = 3.9 and (3 = 4.05 to 
be the systematic error due to cut-off effects. We therefore take r$ = 0.465(6)(14) fm, which is in agreement with ro 
extracted from the value of the pion decay constant, f n , when converting our results at /? = 3.9 and (3 = 4.05 in units 
of r . Using the value of r /a = 5.22(2) and r /a = 6.61(3) at (3 = 3.9 and (3 = 4.05 we find the for the lattice spacings 
c^^.g = 0.089(4) fm and 0^=4.05 = 0.070(3) fm. These values for the lattice spacing extracted using the nucleon 
mass at the physical point are in agreement with those determined from the pion decay constant. This constitutes a 
nice consistency check of our lattice formulation. In what follows we will use the lattice spacing extracted from the 
nucleon mass using 0(p 3 ) heavy baryon chiral perturbation theory to convert the rest of the masses into physical 
units. 



V. CHIRAL EXTRAPOLATION 

Given that the cut-off effects are almost negligible in our simulations, we apply continuum chiral perturbation 
theory to extrapolate lattice results at j3 = 3.9 and (3 — 4.05 to the physical pion mass. In particular, we will use 
SU(2) chiral perturbation thcory(xPT) [47] for two reasons: The first being that our simulations are done for two 
mass-degenerate dynamical quarks and the second because it was shown that SU(3) xPT fails to describe lattice 
data [44]. We would like to stress however that the issue of the applicability of SU(3) %PT is not entirely settled and 
e.g. SU(3) fits to lattice results using staggered fermions were claimed to be produce reasonable fits [48, 49]. 

The leading order SU (2) heavy baryon chiral perturbation (HBxPT) results are given by 

">x°W=4 ) -4 ) < > ( 22 ) 

with two fit-parameters, the baryon mass in the chiral limit 777 ^ and the latter of which gives the leading 
contribution to the dx-tcrm. 

The leading one-loop results for the nucleon and the A in HBxPT were first derived in Ref. [50] and successful fits 
to lattice data on the nucleon and A were discussed in our previous study [4]. A natural generalization of the 0(p 3 ) 
results for the nucleon and A to the rest of the octet and decuplet baryons [51, 52] is given by 

miv(m w ) = 

ra-A("l7r) = 

m^(m n ) = 
7773(777^) = 

for the octet baryons and 

TOa(TO7t) 

"7s* (m„) 



(0) 
N 



,(0) 



,(°) 



4c$ m 2 



4 C W 777^ 



(1) 9 
771 4^' — 4C£ 7774 



,(0) 



4ci 1) m^ 



9 ^ ml 



167T/2 



2.9 



EE 



167T/2 



3 5 1 



167T/2 



(23) 



4 0) 

(0) 



44 1} ml 



4<4^ m 2 

^ L f2 " l TT 1 



25 3aa 3 
— 777 

27 167T/2 - 

9 16tt/ 2 m ^ 

5 3 
3 167T/2 m - 



(24) 



for the decuplet baryons. 
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In addition we consider a cubic term of the following form 

Tnx{fthr) = m ^x ~ ^ C X m \ + C X m Tt (25) 

(2) 

treating c x as an additional fit parameter. 

The next to leading order SU(2) xPT results [47] for the octet are given by 



- H™«, A AS , A) - 0^ I 



m» LO (m«) = m^°(m v ) - ^% ml - T{m^ -A AE , A) - -Mf^ F(m v , A SE ., A) 



167T/2 - 3(4^)2 ^ *' 3(47^)2 

3ffI H 3 2 g |. 5 

■— "— ' 777 — 

167T/2 - (47T/ W ) 5 



mS LO (m„) = m|°(m T ) - to 3 - J 

and for the decuplet baryons: 



2 



to 



(m„-) = m%?{m„) - ^ ff 13 ^ m 3 - 2 [g|. s ^(to^, -A SS ., A ) + gis» -^(m^, -A as »,a)] 

9 16tt/ 2 3(47T/ 7r ) 2 

5 « 2 n 2 

NLO/^ \ _ ^iO/ s .95*5* „,3 #5*5 



3 167T/2 - (47T/ W )2 

< L °(«V) = m^m*) (27) 
with the non analytic function [53] 



A - VA 2 - to 2 + ie \ 3 A 2l fm 2 \ a3i /4A 



^(m, A, A) = (m 2 - A 2 )\/A 2 - m 2 + log V - Am 2 log - A 3 log — (28) 

\A + V A 2 - to 2 + ie y 2 V ^ / V m / 

depending on the threshold parameter A^r = my ^ — and on the scale A of chiral perturbation theory, fixed to 
A = 1 GeV. For A > the real part of the function JF(to, A, A) has the property 

. f -Tim, A, A) m < A 

.Fm,-A,A = ' 2 3/2 (29) 

( —J- (to, A, A) + 27r (to — A J to > A 

which corrects a typo in the sign of the second term in Ref. [8]. In our fits, the nucleon axial charge gA and pion 
decay constant f n are fixed to their experimental values (we use the convention such that / OT = 130.70 MeV). The 
remaining pion-baryon axial coupling constants are taken from SU(3) relations [47]: 

Octet : g A = D + F, 5ss = 2F, g~~=D- F, g AS = 2D 

Decuplet : 3a a = H, gs*s* = §W, g^s-s* = \U (30) 
Transition : g&N =C, = -^C, gz*z = ^jC, g AS » = 

As can be seen, in the octet case, and once gA is fixed, the axial coupling constants depend on the single parameter 
written as a = jj^rp- Its value is poorly known. It can be taken either from the quark model (a = 3/5), from the 
phenomenology of semi-leptonic decays or from hyperon-nucleon scattering. We take 2D = 1.47 or a = 0.58 as given 
in Ref. [47] . The decuplet coupling constants depend on a single parameter for which we again take the value Ti. = 2.2 
from Ref. [47]. This value is not far from that predicted by SU(6) symmetry, H = |<?a = 2.29 used in our previous 
work [4] resulting in the same cubic term for the nucleon and A. For fixing the octet-decuplet transition couplings 
we take the value C = 1.48 from Ref. [53] . 

With the coupling constants fixed in this way, the LO, the one-loop as well as the NLO fits are left with the two 
independent fit parameters and c~x ■ All mass parameters to^ are treated independently unlike what is done in 
Ref. [47] where a universal mass parameter was used for all barons with the same strangeness. 
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A noticeable result of this expansion is the absence of a cubic term in the expression for the A and f2 masses given 
in Eqs. (26) and (27). In the case of fl, it follows from the absence of light valence quarks. However the absence of a 
cubic term in the NLO expression of A, although a consequence of X-PT, is nevertheless a questionable result, since 
it relies on the assumption that m^ <§C Ms — Ma- In the limit A — > the non analytic function (28) becomes 

•F (777*77- , A — > 0, A) = 7T7?7^ , 

which generates a cubic term for the A and slightly modifies the one for S. The corresponding expressions are given 
by 

(°) _aJXI m a _ .9as 3 



m\{m v ) = my - Acy m, - 



167T/I * ' 

2 ffls +ffAs/ 3 „ 
167T/2 



m s (m.) = - fc&'mj - ^ ' ^ , (31) 



in agreement with the results of Eq. (23). 
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FIG. 14: Chiral extrapolation of the nucleon mass in 
units of ro ■ Filled triangles and squares are results at 
(3 = 3.9 on 24 3 x 48 and 32 3 x 64 lattice sizes respectively. 
Open triangles are results at j3 = 4.05. We show chiral 
extrapolations linear in ml as in Eq. (22), to 0(p 3 ) as 
in Eq. (23), fit in the cubic term as in Eq. (25), NLO 
and NNLO in SU(2) chiral perturbation theory as in 
Eqs. (26,32), respectively. We include an error band only 
for the 0(p 3 ) fit for clarity. The physical point shown 
by the asterisk uses the value of ro extracted from f n , 
whereas the cross uses determined from the nucleon 
mass. 
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FIG. 15: Chiral extrapolation of the A mass. The nota- 
tion is the same as that in Fig. 14 but in physical units. 
Here we also include an error band for the cubic fit. 



The expressions for the strange baryon masses to NNLO in xPT given in Rcf. [47] involve in general more unknown 
low energy constants and only if we perform a constrained fit we have enough data to extract them. We found however 
no real advantage in using constrained fits, which generally gave larger x/d.o.f. and did not improve the prediction of 
the mass at the physical point as compared to the unconstrained fits. For the nucleon, A and O masses, unconstrained 
fits can still be performed with four fit parameters [8], namely to^ -*, a and (3 appearing in the expressions of 
NNLO xPT given below. 
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FIG. 16: Chiral extrapolation of the E mass in physical 
units. We show chiral extrapolations linear in mj;, using 
Eq. 25, C>0 3 ), NLO and NNLO in SU(2) chiral perturba- 
tion theory given in Eqs. (23, 26) respectively. The rest 
of the notation is the same as that in Fig. 15. 
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FIG. 18: Chiral extrapolation of the S mass. The notation 
is the same as that in Fig. 16. 
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FIG. 17: Chiral extrapolation of the E* mass. The no- 
tation is the same as that in Fig. 16. 
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FIG. 19: Chiral extrapolation of the H* mass. The nota- 
tion is the same as that in Fig. 16. 
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FIG. 20: Chiral extrapolation of the A mass. The notation 
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and NNLO in SU(2) chiral perturbation theory given in 
Eq. (34). The rest of the notation is the same as that in 
Fig. 14. 



and the real part of J satisfies 

, ( J(m,A,\) to < A 

J(m ^ A ' A) = lj( ra ,A I A)-2,AK-AY /2 m>A (36) 

Using the above Ansatze the chiral extrapolations of lattice results at /3 = 3.9 and j3 = 4.05 given in Tables V-VII 
are performed. In Fig. 14 we show the fits for the nucleon in units of rg. The nucleon mass at the physical point has 
been expressed in units of r$ using the value determined from nucleon mass as well as from the pion decay constant. 
As can be seen these values are consistent. The 0(p 3 ) being the one used to determine the scale passes through the 
physical point. The other curves show the dependence on the chiral Ansatz used. It comes as no surprise that the 
NLO result docs badly for the nucleon underestimating the mass at the physical point whereas the NNLO fits over- 
correct and yield a higher mass. Lattice results at (3 = 3.9 and 4.05 expressed in units of r$ fall on a universal curve 
confirming that finite cut-off effects are small. Therefore we corroborate the conclusion that we can use continuum 
chiral perturbation theory to extrapolate lattice results at j3 = 3.9 and (3 = 4.05. For the chiral extrapolation of the 
other baryons we use the scale determined from the nucleon mass to convert to physical units. 

We show in Figs. 16, 18, 20 the chiral extrapolations for the octet baryon masses and in Figs. 15, 17, 19, 21 the 
corresponding fits for the decuplet baryons given in physical units. We emphasize that the physical point is not 
included in these fits. 

The LO expression describes well the lattice results but leads to extrapolated values inconsistent with the exper- 
imental point. The 0(p 3 ) HBxPT expansion given in Eqs. (23) and (24) with two fit parameters m*- -* and 
provides a good description of lattice data and the results extrapolated to the physical point are in general in agree- 
ment with experiments. The NLO leads to a clear improvement in the case of the A and 2 masses, whereas for the 
rest of the baryons the improvement is marginal. Apart from the preceding remarks, there is no clear advantage in 
using higher order fits, especially NNLO, which even turns out to be numerically unstable for the case of the A and 
Q masses. Therefore our main conclusion is that the 0(p 3 ) HBxPT provides a reasonable description for the nucleon 
and A masses whereas NLO SU{2) for the the strange baryon masses, yielding values at the physical point that are 
in agreement with experiment. 



~Olj?) 64.2(8) 34.7(9) 37.1(8) 9.7(1.1) 62.2(1.1) 38.0(1.7) 17.3(1.3) 6.3(1.3) 

fit with cubic term 33.4(6.9) 33.7 (9.8) 35.6(10.6) 26.3 (7.7) 25.7 (2.5) 24.2 (6.4) 32.4 (13.4) 2.9(1.4) 
NLO 92.5(7) 52.8(8) 43.3(9) 17.2(1.0) 79.5(1.0) 44.1(1.7) 27.9(1.3) 6.3(1.3) 



TABLE XIV: a-term in MeV using the fit parameters determined from 0(p 3 ) xPT, using a cubic fit Eq. (25) and NLO. We 
used the scale from the nucleon mass to convert to physical units. 
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FIG. 22: Comparison of masses for the low lying octet 
baryons. Results from this work are shown by the filled 
(black) triangles for L = 2.1 fm and (blue) squares for 
L = 2.7 fm with a — 0.089 fm and with the open (red) tri- 
angles for L — 2.1 fm and a = 0.070 fm. Results with the 
hybrid action (LHPC) are shown with the (green) aster- 
isks for a = 0.124 fm and results using Nf = 2 + 1 Clover 
fermions (PACS-CS) with the open (orange) circles and 
a = 0.0907 fm. For the nucleon we also show results using 
Nf =2 + 1 asqtad improved staggered fermions (MILC) 
denoted by the filled (light blue) circles. The physical 
masses are shown by the (purple) star. 
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FIG. 23: Comparison of masses for the low lying decuplet 
baryons. The notation is the same as that of Fig. 22. 



We use the relation m£ ~ /x to evaluate the nucleon cr-term by computing m% . Using our 0(p 3 ) fit we find 

crjv = 64.2(9) McV in agreement with the value given in Ref. [4]. Generalizing the relation we can evaluate the 
corresponding tr-terms for the other hadrons. We list in Table XIV the values we obtain using the nucleon mass to 
set the scale. As can be seen, the value extracted depends on the fitting Ansatz. In the most interesting case of the 
nucleon the result of 0(p 3 ) changes by two standard deviations if the coefficient of the cubic term in m n is fitted. In 
the case of the A fitting the cubic term gives the same value as that obtained to 0(p 3 ), compatible to that of the 
nucleon. This is another indication of the argument presented above in favor for the presence of a cubic term in 
of comparable size as that of the nucleon. In fact the main conclusion of this exercise is that allowing the coefficient 
of the cubic term to be determined from the data produces a er-term that for all baryons except the SI is comparable 
within error to the value of nucleon cr-term. Comparing to the results of NLO we can see that for the nucleon this fit 
produces too much curvature as already observed for instance in Ref. [8] . For the other baryons a reasonable value is 
obtained depending on the quality of the fits as seen in the figures. 
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VI. COMPARISON WITH OTHER LATTICE RESULTS AND WITH EXPERIMENT 

In this section we show a comparison of recent lattice results on the baryon masses from various collaborations in 
Figs. 22 and 23. For our results we use the lattice spacing determined from the nuclcon mass to convert physical units. 
Results from the other collaborations are converted to physical units using the lattice spacing that they provide. The 
level of agreement of lattice QCD results using a variety of fermion discretization schemes before taking the continuum 
limit or other lattice artifacts into account is quite impressive. Small discrepancies seen mainly in the decuplet masses 
can be attributed to lattice artifacts. In particular results using asqtad improved staggered fcrmions may suffer the 
most from discretization errors. The MILC collaboration has simulations on finer lattices and an update on the 
masses is expected in the the near future. We note that results very close to the physical point obtained using Clover 
fermions from the PACS-CS collaboration [6] may have large finite volume effects due to the fact that m^L < 3.5 in 
this simulation. 

Finally we show our continuum results in Fig. 24. We take the continuum limit using results at /3 = 3.9 and 
(3 = 4.05 after interpolating at a given value of r^m^. The continuum values used are collected in Table XIII and 
VIX for the octet and decuplet respectively. Residual cut-off effects that may result from using Eq. (10) to estimate 
fx s at (3 = 4.05 are not included in the systematic errors. For the nucleon and the A we use 0(p 3 ) to extrapolate to 
the physical point as done in our previous work [4]. For the strangeness non-zero baryons we use NLO SU(2) HBxPT 
to extrapolate to the physical point. In the statistical error we have added the error arising from the uncertainty in 
?o , i.e. the difference in the resulting masses when we use r$ = 0.471 fm and r$ ± 0.021 fm to set the scale. As can 
be seen, our results compare well with experiment within the estimated uncertainties. 

Baryon Spectrum 



I" 

i_ 



Exp. results 

Width 

* Input 

• Lattice results 



12 3 

Strangeness 

FIG. 24: The octet and decuplet spectrum. Data are shown in the continuum and at the physical point. For the nucleon and 
A using 0(p 3 ) and for the rest using NLO SU(2) HBxPT for the chiral extrapolation. 



VII. SUMMARY AND CONCLUSIONS 

The focus of this work is the computation of the low-lying baryon masses using twisted mass fermions at maximal 
twist. It is in line with ongoing efforts by other lattice collaborations worldwide to predict the baryon spectrum from 
first principles. Comparison of lattice results with experiment is regarded as an important benchmark for lattice QCD 
and justifies the use of different lattice actions, each with different systematic errors. For example, the twisted mass 
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action with only one parameter to tune, provides automatic 0(a) improvement. However it restores isospin symmetry 
only in the continuum limit. We have examined the consequences of isospin breaking on the baryon masses and found 
them to be either small or compatible with zero. On our finer lattice at ft = 4.05 the maximal isospin violation 
is obtained in the octet only in the case of the Cascade where we find that m^o — m s - ~ 50 MeV. Finite volume 
corrections are estimated in the case of the nucleon and found to be smaller than statistical errors. The continuum 
extrapolation using results at ft = 3.9 and ft = 4.05 are verified using a finer lattice at ft = 4.2 in the case of the 
nucleon supporting the analysis carried out. Therefore this study shows that both cut-off effects and finite volume 
corrections are small and continuum results can be extracted using lattice data at ft = 3.9 and 4.05. 

An investigation of the Gell-Mann-Okubo relations has been carried out. For the octet baryon masses we find that 
these relations are satisfied at all pion masses even at a non-vanishing lattice spacing. For the decuplet baryon masses 
we sec deviations and it will be interesting to study these relations at finer values of the lattice spacing and smaller 
quark masses. 

Comparison of the results at given lattice spacings with those of other collaborations reveals consistency among 
groups using improved actions with lattice spacing being smaller than 0.1 fm. This is a non-trivial consistency check of 
the several lattice formulations directly on lattice data without without the necessity of any extrapolations. This level 
of agreement among different lattice actions, is a welcome outcome of the collective effort of several collaborations. 
The final continuum results at the physical limit shown in Fig. 24 are in excellent agreement with experiment. The 
largest uncertainty in the final value comes from systematic errors in setting the scale, which are an order of magnitude 
larger than statistical errors. 

Besides the masses we have extracted from the various fits the er-tcrm. To 0(p 3 ) x?T we find a value of 64(1) MeV 
for the nucleon er-tcrm. Allowing the coefficient of the cubic term in to be determined from the data yields a 
smaller value of 39(12) MeV albeit with a large statistical error. Fitting with the latter Ansatz produces for all 
baryons expect the f2 a value of the cr-term compatible with that of the nucleon. Clearly this is a result that requires 
further study. In particular results at smaller pion masses will help to better determine the curvature of fits. 

The next step for the ETM collaboration is to perform the analysis using a dynamical strange quark. Within the 
twisted mass formalism this is accomplished by simulating a non-degenerate doublet. Such Np = 2 + 1 + 1 simulations 
are already available at two values of the lattice spacing [54] comparable to the ones studied in this work. This future 
study will provide a nice comparison to the present work and enable us to gauge unquenching effects in the strange 
quark sector. 
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